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Soife  ,'i.i a? >ut:\ t:  Lnv» at  f'roc aduret.  for 
Ss  ii  ' ? tin;,;  a Cumulative  t'istvibat.ion  Function 

1.  Suna-ary.  tome  invariant,  procedures,  which  are  essential!}  step  functions, 
are  considered  as  estimators  of  the  cumulative  distribution  function  of  a uni  - 
variate  random  variable  on  which  a finite  fixed  number  of  observations  are 
given,  for  various  loss  functions.  Two  principal  classes  of  loss  functions  are 


considered  and  it  is  shown  that  for  a special  loss  function  in  one  class  the  op- 
timum orocedure  is  the  usual  sample  r-nmiT j> tu-so  function. 

2.  Introduction.  Suppose  that  a sample  Xj,  ^2’ ' ‘ '^n  a one-dimensional 

chance  variable  X is  given.  In  a recent  paper,  Birnbaum  B has  discussed 
various  techniques  for  deciding  whether  X has  a completely  specified  continu- 
ous cumulative  distribution  function  (c.d.f.),  H(x)=sP(X^x) . In  this  paper 
is  discussed  an  allied  problem,  viz.,  that  if  F(x)  sb?(X^s)  is  the  unknown  con- 
tinuous c.d.f.  of  X and  if  F(x)  be  an  estimate  of  F(x)  based  on  the  sample 


X,,X_, . . .X  , what  would  be  the  best  estimate  F when  certain  forms  of  the  loss 
12*  n' 

function  are  given. 

Consider  the  loss  function 


.-CO 

(1)  L(F,F)  = f |f(x)~F(x)  [ dx 

— 03 


where  r is  an  integer  JSrl.  It  is  almost  obvious  that  the  only  invariant  pro- 
cedures for  estimating  F under  the  group  of  all  cne-to-one  monotone  transforma- 
tions of  the  real  numbers  onto  themselves  which  leave  the  sample  values 
Xj_  (iP=l,2, . . .n)  invariant  are  those  which  esMmata  F(x)  by  a step  function 

(2)  F(x)  = constant,  say  c^  for  X^^  x<X^f1^ 

« 0 

where  X(L)<X(2)<  . . . <X(n) 

are  the  ordered,  observations  and  and  X^n+x^ 


denote  -i*»end  ±co  respectively. 


Gsing  this  entira:  ,e  F,  we  got 


L(F,F) 

j=-0 


vOt1) 


|f(x)-c. 


T dF(x) 


(3) 
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where  X-,  <Xn  is  an  order'd  sample  of  sixe  n from  this  rectan/ular 

distribution  over  (0,1) ,Xq  and  Xn^i  denote  0 nnct  1 respectively,  ond  the  symbol 
E denotes  that  the  expectation  is  taken  with  respect  to  the  rectam-ulor  distri- 
bution over  (0,1).  In  the  rest  of  this  paper,  we  shall  use  consistently  the 
letter  E to  denote  the  fact  that  the  expectation  is  to  be  taken  with  respect  to 
the  rectangular  distribution  over  (0,1). 

The  same  argument  applies  when  the  loss  function  is  of  the  form 

(5)  L(r,h=  dF(x) 

J F(x)  [l-F(x)J 

—co 


3 


A 

and  in  this  case  by  taking  V as  in  (.?;  we  obtain 


(6) 


R=R 


±_  f*j+i  h^dl 

T*0  \ *(!-*) 


dx 


where  X-j;  j = 0,l, , . .n+1  are  the  same  at  in  (4). 

It  is  obvious  that  since  risk  R is  constant,  a minima*  procedure  will  be 

to  choose  Cj; j— 0,1. . .n,  such  that  R is  minimum.  We  consider  in  this  paper  the 

values  of  Cj  when  the  loss  function  is  of  the  form  (1)  for  all  integers  r S^l 

si  1 and  when  r is 

and  when  the  loss  function  is  of  the  form  (5)  for  rA*n  even  integer  ^2  The 
case  when  r is  odd  in  (5)  seems  to  be  rather  complicated. 


3.  T he  lose  Amotion  L(F,$)=  j [f(x)  - F(x)]2  dF(j) 


In  this  case 


(7.) 


a.  rl 

R=SS^Jx  (X'CjrdX 
•• 

=1/3  Efe-X3>  - 3cJUH1-X3)  t^(xJ,rxp] 


Since  the  distribution  of  the  j-th  order  atatistic  Xj  in  a sample  of  sioe 
n from  the  rectangular  distribution  over  (0,1)  is  a beta  distribution  with  prob- 
ability density 

(8)  p(y)  - -^7^-777  yJ"_1  (i-y)°-J  , osiyil  , 


It  is  easily  scan  that  for  any  positive  integer  r. 


(9) 


E(XT)  — ^ • n ’ 1 

J (n*l)(n*2) (n>r) 


c io)  e(x?,  ~iV)  .Utiini). 

J r"  |n4l) ! iif2) — (n+r) 


•^1 


TTT 


f c r r 1 . 


It  will  be  u^e^ul  to  remark  that  (10)  holds  for  all  j;  .5*0.  I,. 
Substituting  in  R we  obtain  after  some  simplification, 


. ,n. 


(11) 


R ss  1/1  - 


- 1 


« i (j^  +iti  f:  **■ 

J=U  j»0 


(n-»-l)  (n+2) 


1~  (CfM)2 

n+1  4-  1 


6(n+2)  ' n+1  2L-.  v~j  n+2 

j*U 


We  see  thus  that  R in  minimised  by  choosing 

(12)  cj=  5 j"  0,1,... n. 

and  hence  the  minimax  invariant  procedure  is  to  estimate  F(x)  by 

(13)  ^x)sS"^  ? XJil<X^i  ? Js0,l,...,n, 


where  (Xj  ,X,,,, , . X^)  is  the  ordered  Beraple  and  X^  and  X^j  stand  for  -00  and  fcz> 
respectively. 

The  minimum  risk  corresponding  to  this  procedure  is  seen  to  be  l'/’6(nf2). 

It  is  of  some  interest  to  note  that  the  risk  corresponding  to  the  usu  c pro- 
cedure of  taking  0^=1  j 'tx  is  given  by  1 ^6n 

00 

]F(x)  - f(x)j  dF(jv)  . 


4-  The  loss  function  L(F,F) 


- 


For  tnis  case 


(U) 


3 


I 


vnere 


= 


r1  ; 41 

F,l  * jx-Cijdx 
j 


* H 


(15) 

Now 

(16)  E 


f VCj)  / XrCj)]=  I (t-'j)2  p(y)d>-  -f  j (y-Cj)2  .jlyjdy 

Jc.  -'0 


where  p(y)  is  Riven  by  (8),  ami  sinilarly  we  can  Ret  E [”(X  -c  )/x  -c.O  = 

**  jfl  j ' jfl  JJJ 

Substituting  in  (15)  we  easily  obtain 


/o  , i (Hi 


r1 

(a-j)  (y-c,)V  (1  -y)n_j“‘L 


Jc 


dy 


(17) 


- (o-J)  I (y-c,)2  yj  (l-y)n""^_1dy 

J0  J 

■ - [ (y-Cj)2  yJ"1{l-y)n_,i  dy 

fCi 

c-1 


r 

t J (y-c,)2  r!”1(i-y)n“'3  dy. 
Jo 


This  eventually  leads  to 


. o a r n-j  Q Ji-kf2  I 

(IB)  rj  = (j)  - .jXHX.-HD  tag  <-«V;b 

for  j S 0,1,2, . . ,n. 

Since  R«  ^L.  ^ , and  fro~i  (15)  we  see  that  for  each  j,£j  is  positive  and 

iepeads  only  on  it  is  obvious  that  to  minimize  R,  it,  is  necessary  and  suffi- 

/? 

cifetil  to  minimize  each  <*»  ^ separately.  we  have 


6 


d Si,  r n-i  • Or  f 

(19>  gsj  = 9 <j'  ♦ gj  <-«*<V>  -tsa-J 


cj 


= 2 (“)l  | *4(1-2)*- J *3-  SL(j+l,n-j+l) 

L: 


and 

(20) 


= 2(j>  °i  ‘W1 

J 

95^ 


Setting  9c,  = ° and  solving  wa  obtain  cj  as  the  median  of  the  Beta  distribution 

j 


with  density 

(21)  g(»)  - 


B(  2 rl»o— jf 1) 


(1-3  n'j  , 0 1 , 


for  j*0,  1,2, . . . ,n. 

Since  (20)  shows  thao~— ^>0  for  0 ^ c , <.1,  it  follows  that  this  solution  for 

cj  ' 

cj  in  fact  minimizes  dj  for  . . . ,n,  and  herce  minimizes  H.  The  minlmax  in- 

variant procedure  is  thus  t)  estimate  F(x)  by 

F(x)  - c . X r£  x < X 

i i Jil 

j^O,  1, . . ,n, 

where  (X^Xg, • ■ . ,XQ)  .La  the  ordered  sample,  Xq  end  X ^ stand  for  - co  and  f co 

respectively,  and  c,  ( j 0,1, . . . ,n)  is  the  medisn  of  the  Beta  distribution  with 

density  (21).  An  alternative  way  of  obtaining  this  result  is  piven  in  section  9. 
It  is  rather  interesting  to  note  that  for  the  loss  function  discussed  in  the  last 
Section,  cj  was  obtained  as  the  mean  of  the  same  Bata  distribution. 

The  actual  computation  of  the  values  of  cj  (js0,l,...a)  can  be  easily 

carried  out,  for  a given  n,  with  the  help  of  the  tables  of  the  incomplete  Beta- 
function  t 3 I • In  the  notation  of  the  tables 


7 


f X P-A  (l-xj'^dx 

, . . , Jo 

(22)  QJ  = 

Z1 

\ xP-i-d-xJ^-idx 
) 0 


Thus  we  have  to  find  the  value  of  x such  that 


(23)  Iy  (jfl,n-j*l)*  3 

Usinf  the  relation 

(24)  Tx  ( p,a) ~ 1-It_x  (hjp)  • 


it  is  seen  that 


(25) 


C_  , = 


'n-j 


1-e 


and  thus  only  spoilt  half  t.ne  total  number  o'1  r values  have  to  ! r actual  1;  o’! 
tained  *Yon  the  tables  The  values  of  ct  (j  s 0,1,  ,n)  for  n *0*1,2.,  P-2 

to  tw0  •iecirwl  .-laces  have  been  co:^  ut.ed  and  labui'  ter*  below. 


.Tiible.,! 

Values  of  cj  ( j =0,1, . . . ,n) 


for  n«  1,2, 12 


c0 

■ 

c2 

c3 

" 

c4 

c5 

r 

c6 

7 

c7 

p3 

e9 

7~T 

' 10 

~T 

cn  i 

el? 

1 

.29 

.71 

| 

n 

.21 

• 50 

.79 

3 

.16 

.39 

.61 

.84 

4 

13 

• 31 

.50 

.69 

.87 

5 

.11 

.26 

42 

• 53 

.74 

89 

6 

• 09 

• 23 

36 

50 

.64 

.77 

• 91 

7 

.08 

.20 

.32 

.44 

• 56 

.68 

.80 

• 92 

8 

.07 

.13 

.23 

• 39 

• 50 

.61 

.72 

.82 

• 93 

9 

.07 

.16 

.26 

.35 

■ 45 

• 55 

•65 

-74 

.84 

• 93 

10 

.06 

.15 

.23 

• 32 

.4-1 

• 50 

.59 

.68 

.77 

• 85 

• 94 

11 

.06 

.14 

.22 

-an 

.38 

.46 

• 54 

.62 

.70 

.78 

.86 

• 94 

12 

L!i 

_±1 

.20 

.27 

• 35 

.42 

• 50 

• 58 

.65 

.73 

.80 

.87 

.95 

...  . _ 

c> 


5 The  lo^?  fun 


sao 

ction  L(F, F)  - J |f(x)  - F(x)}  2/f'(x)  ji-F(x)  j dF(x) . 


For  this  loss  function,  ns  mentioned  before,  we  pet 


J=0 


r-X  j t t ...  - >2 

JtA 


J X, 


x(l-x) 


dx 


(26) 


_ (h  (*-c0): 


“•f 

■J  0 


x(l-x) 


n-1 

dX  S*  E 1 ~^XJ+i"V  * cj2(logXjti-loexj^ 


-d-cj 


)2£log  (1-Xj+1)  - lop  (1-Xjjjj+  e( 


Crv' 


x(l-x) 


x2 


For  finite  risk  the  integrals  in  the  first  and  the  last  terms  of  the  above 
expression  must  be  finite.  The  necessary  and  sufficient  condition  for  this  is 
that  Cq  — 0 and  cn~-l«  Our  set  Cj  (j  = 0,l,...,n)  must  then  be  such  that 


Cq  =*  0 and  cQ  * n . 


With  the  convention  that  0 :c  co  = 0,  we  cfin,  therefore,  write  (26)  in  the 


form 


(27)  R*  21  eFI<X  -X  ) i c.2(iog  X -log  X.)-(l-c.)2(log(l-X,  )-log(i-X 

J=o  | JT-t  J .!  Jr1  J J C_  Ji-1  J 

The  probability  density  of  X is  given  by  (8),  from  which  we  obtain 

J 

c1 

(28)  E(log  X)  » j (?)  7J_1  (l-y'S  i log  y dy 

J Jo 

and 


(29'  Sfloe:  (1-X,)»  j (?)  f yj  3 

4 J \ 

Jo 


(l-y)”  3 log  (l-y)dj 


In  order  to  evaluate  (28)  and  (29)  we  make  use  of  the  following  lemmas: 


S' 


(30) 


1 ( l-y)n^ 


log  y dy--:  (XUBfUJ.IL  fy(j)-)Ka+l)]where  y'(k)=p'  (k)/p(k) 
P(n-tl)  L -3 


O*  f y^’^d- 


y)n~j  dy.  The  left  hand  side  of  (30)  is  f 1 (c* ) 


Proof.  Let  V*,-  f 

J0 

evaluated  at  CV-~j  as  can  he  seer.  by  differentiating  under  the  integral  sign. 

But  f (oS)  “ P (°^  ) P(n-j|l) -//'1(0<+  n~jfl),  and  the  desired  result  is  ob- 
tained by  evaluating  the  logarithmic  derivative  of  f(c?<)  at  e(  & 

Lemma  5.2; 

(31)  C yj”1  (l-y)n~j  leg  (j.-y)dy z^CLAllL jjt'(n--jil)  - (n+l)J 

v 0 

where  J^(k)  — p*(k)/p(k). 

Proof.  Exactly  as  in  Lemma  51,  or  easily  obtained  from  it  by  a change  of  variables. 
Utiliz.i  ig  Lemmas  5-1  and  5 • 2>  we  obtain 


(32)  B(lcg  Xj)  * f{i)  -^(a+i;- 

and 

(33)  E log  (1~X.)~  ^(n-jfl)  -J^(ihI), 

J 

where  /'(k)  ==  fV  (k)  / p (k) . 


Further,  since  p(kfl)  ax  kp(k)  , p*(kfl)«  P(k)  ikp*(ic)  , 

we  r.se  that  )f'  (kfl)  — p*  (k-rl)/p(k+l)~  l/k  -5- ^(k) , 
and  hence  the  function  ^ satisfies  the  difference  equation 

(34)  y*  (k+1)  -J*(k)  — l/k. 

From  (32),  (33)  and  (34)  we  get 

(35)  E(log  X.fl  - log  Xj)—  l/j  , for  j O 
an'’ 

e[ log  (1-X.  ) - log  (1-X  ) -l/(n-«)  , for  J^n. 


(36) 


Subr.lif  .titi-  from  (20),  (35)  . and  (36)  .n  (27)  we  obtain  at  once 

, , ? .SrJ-  f ' -l  o i .oli  n-i  » . . 

(37)  R -r-;r'  -!•  - ---  -r  •?  cf-/-  (1-e.r  L.±  “ (e- -i;' 

fi(nrl,  j*i  n+i  .i  J n-j  ,1  J n Jt.n-J)  2 n' 


It  is  seen  from  (37)  that  R is  minimized  by  choosing 


(38) 


i j =5  j for  j*l,2, 
Since  Cq  •“  0 and  cQ  1, 


. (n-l). 

this  expression  for  holds  good  also  for  j — 0, 

J 


and  j — a.  Tons  the  nrinimax  invariant  estimate  F for  the  loss  function  in  this 
Section  turn/,  out  to  be  the  usual  sample  cmul.atJ.ve  function 

when  X^  **  , j=0,l, . . . ,n. 


(39) 


F(x)< 


c.^j/n, 


where  X^<  X2 < ....  <XTJ  is  an  ordered  sample  from  the  c.d  f.  F,  Xq  and  Xn+1 
standing  I or  - «>  and  fco  respectively.  The  actual  value  of  the  risk  corres- 
ponding t.D  this  estimate  is  l^n. 


6.  The  loss  function  L(F,^)a J |f(x)  - F(x)| /F(x)£i-F(xj]  dF(x) . 
~<oo 


In  this  case  we  obtain 

(40)  R E ]>_  \ jx-c^  j /x(l-x)dx  - 'y  fit-  > 

j-0JT  j J=0  J 

■ j 


where 

/'X  . ., 

(41)  j x-c j j /x(l-x)  dx 

“/X, 


As  in  the  last  Section,  it  will  be  seen  that  for  finite  risk  the  necessary 

and  sufficient  condition  is  that  c,,— 0 and  c s=-l . For  j*£o.  n.  we  obtain 

0 n / 


l.i 


(42)  & = 2 


’} K Vloe  *i  l-«il lof‘  vl0i-  V I i(1-c})l:l K>cU-=j)-iog< 


-ti-Oj)  j log'1-0.)  - log(l- 


:1  xJ>j] 


The  distribution  of  Xj  hap  probability  density  p(y)  g iven  by  (8)  and  tb 
distribution  of  X.  ^ has  the  probability  density 

(4i)  q(y)*  yj  (l-y)1 *^-1  ,0^y<l. 

-“vj  t ) 

Using  (8)  and  (43)  we  can  express  in  the  form 

J 


(44)  Si  = (J) 


where 


C j)(c^y)  dy  ^(0^ ,y) 

Jo  c -V 


(45)  ^(Cj,y)=>  ^cjlog  Cj+(1-Cj  )log(l-Cj ) - c^log  y^l-c^logCl-yjjy^-3  (1 


Straightforward  integration  leads  to 


(46)  ^p^(cj*y)dy  s yJ(l-y)n_^|cj(log  Cj-log  y)  \ (l-c^) (log(l-Cj)-log(l- 

-j- ^(c  ,-y)y^-^  (l-y)n“j_1  dy  constant 

0 

which  enables  us  to  obtain  as 

7cj  r1 

\ (cry)yj_1  (:--y)n"'3_x  dy- l (cj-y)yJ  "(l-y)n~J  1' 


(4?)  a*(a) 

J J 


for  j - 1,2, . . . , n-1. 

Since  $£  and  are  fixed,  and  each  (g  .is  positive  and  depends  only 
On  j 

Q 

to  minimize  R it  is  clearly  necessary  and  sufficient  to  minimise  ^ 

i Ua 

J * 


see  that 


<»>  Ef^  (2Tm-'J)2s41'V.1-  -I.1 

k«0  J • J 

for  j=  0,1,2, ...  ,n. 

Substituting  from  (10)  in  (52)  we  obtain 


T 2s  + 1 (2sfl)(_c  t2sfl-k 

n+I  1 k n cy 


c”  + 


(nfl)  . • (nfk) 


C3) 


1 

nfl 


,»<;a  . 

C1  T 


2sfl 

2z: 

k*2 


(2d  ) (-c,)2s|1-k  Hi-1/1  .:.-y±k--?.Ll 

j (afl) . . . (nfk) 


For  conciseness  we  introduce  the  following  notation  somewhat  similar  to  the  bi- 
nomial and  distinguished  from  it  by  ar.  asterisk: 

(54)  (t  - $ )<*=  ts  r:  (-i)k  (2) 

for  fixed  real  a and  b and  a positive  integer  q. 

It  is  easily  verified  that  for  any  positive  integer  r, 

j* 

(55)  ~dt?  ^ ~ bfl  * q (q-1) (q-rfl)  (t-iili  /q  r when  r*q 

a*  0 when  r>q. 

Using  this  notation  we  can  write 


(56)  * nir  (cj-  -Jji) 


2s* 


We  have  to  choose  cj  so  as  to  minimize  R.  Obviously  minimizing  R is  equiv- 
alent to  minimizing  ^ separately  for  each  j . We  obtain 

ten)  55L-J 2S  jfl,(23-l)* 

n+l  (cj  n+2'  * 


and 


(58)  &L£l  --  , _ Jii.  (2s-Z)  « 

l5J  ?;?  mi  t =i2> 

j 

5 rtf1  , >2s 

= E \ (x-cw  dxXS,  it  is  cl' 

J Sv  f1 


Since 


ear  that 


J 


(59)  = 2s^~l) 


(x-G.)2^  2dx^>0. 


I-et  f(c  )*  ■2“**-  . It  is  easily  seen  that  f(0)  Is  negative  and  f(l)  is  positive, 

J 

and  since  f 1 (c  )2>0  for  all  real  c , f(c  ) is  a strictly  increasing  function  of 

i i i 

c . Hence  f(c  )=-0  for  one  and  only  one  real  value  of  c , ard  this  c neces- 

j j s j j 

sarily  lies  between  0 and  1.  Thus  we  find  that  f and  henro  R,  is  mj.nirnized  by 

dg.  ^ 

setting  g— s 0 and  solving  for  e . the  resulting  equation 
j J 


<ei~iw)<r~1,*=0 


This  equation  has  one  and  only  one  real  root  which  lies  between  0 and  1.  The  min- 
imax invariant  procedure  for  the  loss  function  of  this  Section  is  thus  to  estimate 


F(x)  b3 


F(x)  = Cj  ; X.  ^.xCXj+1 


j«0,l,.  -n, 

where  X.  ; j =0,  1, . n-j-1,  have  been  defined  earlier  and  c.  is  the  real  root  of 
J *j 

(60) . It  can  further  he  seen  from  (60)  that  the  equation  remains  unchanged  if  we 

replace  j by  a-j  and  c.  by  1-c  . Hence  c a*  1-c  . and  we  see  that  in  practice 

J j n-j  j’ 

the  number  of  equations  to  be  solved  is  about  half  the  sample  size. 

It  may  be  roticed  that  for  r — 2,  the  equation  (60)  reduces  to  a l.aear  equation 


which  has  the  uniaue  solution  a,  - as  obtained  earlier  in  (12). 

^ n+ d. 


CO 

8.  The  loss  function  L(F,F)=*  ( Jf(x)  r/F( x)£l-F(x)JdF(x)  where  r is  any 

' ~a> 


tiositiv-:.'  even  integer. 


Lot  c» 2s,  then 


(62) 

where 

(63) 


n (* J+l  (x-cj 

R--E  21  1 — 

j^O  l x(l-x 

V>X_. 


5q 


dx  ~ 


x(l-x)  j=0 


, . / \ 2 s 

( if!  (x"c1) 

E \ — ■ — ” dx  . 

J 1 x(l-x) 

-J  X 

j 


Since  Xg=0  and  X s»l,  it  is  clear  that  in  order  to  obtain  finite  ri3k  it  is 
necessary  and  sufficient  that  Cg~  0 and  cn“l-  For  j^O,  n,  vs  can  write 


/?  , E 


(64) 


where 


2 s-2  1 


isr  hti  % <*5,1  - *5a>  * <=rci°B  V lop  V 

-(i-c.?S  flog  (1-X.  ) _ lop  (i-X  a 

J <-  Ji1  Jjj 


(65)  Iff  (^X-O1  ; h=0,l,2,...,  (2s-2). 

i-0  J 


Substituting  from  (10),  (35)  and  (36),  we  get 
(66)  ( jfh) In!  . 1 2s  . 1 n v2s 

<66)  J & I^hTijTJi  n T 7 J T^(1-:J) 

and  substituting  from  (65),  we  can  write 


(67)  © 


_ 0^-0  . . 9o_2-h 

47  __  n!_  v*  (Jth?!  < -- 


l i+1, 2s.  i,  12s.  1 


, - - > (-D  (7)  e + tc  t “(i-c  ) 


2s 


1=0 


This  ic  2st-h  degree  polynomial  in  c . - Collecting  the  cceffici«'.’;  s ct*  like 


nowers  of  c . 

J 


we  obtain 


16 


_ 2s-2 

//r»^  £ -•  11  2s  2s  2s-l.  \ ^ 

<63)  V 7*537  «j  "IZJ'j  + £o-Vj  ’ 


where 


(69) 


= <-«k11dfl  X 

L ha° 


2a-2-k  “1 

£i  > (Jf*0»  __1_ 

,r  ' (a+h-fl)!  n-j 


for  k— 0,1,2, ...  , 2s-2. 

Tc  simplify  (68)  further,  we  state  and  prove  the  following  lemma: 
T.emmfl  8.1.  If  j and  n are  positive  integers  and  j<n,  then 


n!_  9 (jt-h)!  _ 1 r _ 91.1  jwc 

(n+h+1>!  n~-5  L ot  Jin+rtJ  ’ 


Proof.  The  left,  hand  side 


(j+h)!  (n-j)1. 


J h*0  iwn  n ) ' 


- W ^fh  (1-x)  n_j 

h»0  ) 

-'n 


dx 


= (“)\  {J-^1  ) (1-x) 

-'o 

=.the  right  hand  side,  after  simplification. 
Substituting  in  (69)  from  the  lemma  when  qa«2s-2~k,  we  obtain 


2s-l-k 

(70)  gfe  =■  (-l)k-~  (2S)  FT  ^ k*0,  1,  2,...,  2s-2, 

<*,=1  T 


and  substituting  new  in  (68)  we  obtain 


r 

i 


2s--l->  - 


(7i>  S^M°r  + S * XI  <v*,a" 

in  the  notation  introduced  in  (54). 

Nov  with  the  same  reasoning  as  in  the  last  section  it  will  be  seen  that  £ 

s& , 1 

and  hence  R is  minimized  by  setting  ^ — =0  and  solving  for  c,  the  resulting  e- 

9°j  3 

quation 

(72)  (Cj-  j/n)(r“l)#*0. 

This  equation  has  one  and  only  one  real  root  which  lies  between  0 and  1.  Since 
for  j “0,  (72)  reduces  to  c~r~-=:0  giving  cn  — U as  the  only  real  root*  and  for 
j=sn,  it  reduces  to  (c^-l/'*-^  0,  giving  cn  sz  1 as  the  only  real  root  it 

follows  that  we  can  say  that  the  minimax  invariant  procedure  for  the  loss  func- 
tion of  this  Section  is  to  estimate  F(x)  by 
/\ 

F (x)*  c j j X<s<X.  . 

J J jf-L 

j*«0,  1, . . . , n, 

where  ; j-=0,  1,...,  n-fl  have  been  defined  earlier  and  Cj  is  the  real  root  of 
(72) . Again  the  number  of  equations  to  be  sol  *ed  in  practice  will  be  about  half 
the  sample  size  since  it  can  be  easily  seen  that  (72)  remains  unchanged  by  re- 
placing j by  n-j  and  Cj  by  1-Cj  , so  that  cn_j  =■  1-c^ 

For  r=2,  the  equations  (72)  reduces  to  c . - j/nnO  giving  c.sj/n  for  all  j. 

• J J 


9.  The  loss 


roo 

|F(x)-F(x)  I rdF(x) . where  r is  any  positive  integer. 


J-oo 


In  this  case 


(73)  R = F ^ f ^ ix-cj1 


where 


L 


•*& 


(74)  S.  = -L-  3 

J rj-1 

Using  (8)  we  obtain 


(xj-ri  ~cJpJtr  cj  r - (};rc^lxr 


r!  j' 


(75) 


(Xj  cj>  Jxi_cT  (y-cj)r*1  y'*  1 U-y)11-^ 

~-j 

1 (i-y)n“'i  dyj. 


- ' J(y-yjr*1  y^1 


and  similarly, 


(76) 


:jo 


■1 


E'<V  -cj>  V 


=(n-j)  (7) 


rr1 


1 


-cj 


(y-c  )rt1  y^ (l-y)n  ^ x dy 

J 

1 

(cj-y)r+1  y^(i-y)a_j_1cijr  j • 

J J 


From  (75)  and  (76)  ve  obtain 

[f1 

Uc3  * 

(77)  J 


jf  f / ..  _ /n  / .t  \ Jt_. 

-j  r^l  vj'j  ^ W-j;uy 

c) 


f (-!)  \ (y-c.)r  y^  1 (l-y)n  ” tny~j)dy 

J 


JO 


Again  it  is  obvious  that  to  minimize  R is  equivalent  to  minimizing  !&  for 

J 

each  j . Further  we  see  that  the  conditions  for  differentiation  with  respect  to 
under  the  integral  sign  in  (77)  are  satisfied,  and  we  obtain 


r1 

\ , vr  j-1 

\ (y-c  > r1 


(78) 


LA. 


(l-y)n~^  1(ny-j)dy 


i urf  (y~cj)r  yj_1  (i-y)a"^~J1’(ny-4)dy  \ 
-o  .j 


and  for  r ^2, 


• «r.  , 


&T?  ir 


ip  I ty  T l r3  1 


y-e, 


L»#3 


+ (y-« ,)  r j "1  i l-y } * («;,'• 

2 

'o 


rjfi  , t 

if  <**-!,>  r|  iK^jF  •?  -in>ir 


ItS  Hi 

,f«s*  al*  -£«rir  ter; 


.'Tl'l 


«Vi*  J 

* j 

3 K<XC&  4|  if  Vf*  i **f‘%  CV  ' f«*~  i?n  Y »*  4jnr*  j/**.  Jp^; 


% v?  r**  s 


^ fcr±,b* 


ff  CW.  tt  * ^ 

Jl*  tj-c;  * 


JL  .«,  _ n^srifV'j 

in^i»  ■ nfsu* 


«wl  r*r  *1*.  *>-'1  *y  7,r*  f ■ 1 * la-'nrrS*:  ? 


of  ei , « name*-  ut$  vstlo-  »»  for  <?o*  •jsu  «tw*  r-«X  *■*’•>•*  f \> 


IMp 


fO 


» Tirirr’-rl^f-  i!<r  Nh*  tea  •**■*  "•■a  fir;?  &rnee 

J $& 

IE  Is  «4  fcy  swi-tlr’.  g£*  *£  i4M  far  «j  to*  r "W-ittf  «•. nation 

J 


Cy-c^V^ti-yl0'44 


O 

^y+(-l)rj  | 


'0 


-V J 


r J-i„ 

“ U>tf) 


(ny-i)ay-O 


Thu  a tUr  p&Mm  Ytotuctn  to  tir*&  of  snXXar:  Ifea  «H»w  for  j=r  (%!,,,  , ,n. 

The  ^eoerai  /elution  of  (80)  eivi«f  c,  w*llcitty  lift  terr.*  of*  i*  "#  ^ ^ ilot 

aefia  to  he  0*sible.  We  shf-ll,  however,  sfny  il>  too  e.«n*!oa  *>  th»fc  it 
no^  he  toe  XlfrtcoU  W Obtain  the  solution  in  way  riven  or  so-  It  «*n,  farnimrj 

he  provec  -from  (80>  thsfc  c ■&  1ms  bo  «n*t  lute  ousaher  of  a.  nations  to  b*  solved 

n~3  3 


In  pMettce  Will  ha  shout  half  the  wrJt.de  hi?* 


We.  "fin  fT'itt?  (fu) 
n 


(31)  { '>“^j)V_I(l-:-)n‘J’1(ny-«1)dy^^l-(-l)r^  J (y-e*)V 

~J 0 '-Jo 


(l~y}n0‘A(nJ-:j)dy 


It  would  he  -fen  that  the  rir’nt  hr  nd  Hde  of  till *3  e u.-.tion  would  he  of  H 
oov.ar.;e  only  wrier  r is  odd,  for  when  r is  even,  if,  ?*e  mc°f  tc  »,t.r©» 
irse  left*  h'w.d  Hue  o*’  <?  irif  sen  (til)  ! -.*  f • »r' 1 •" 


(Sk 


fe>  “ 


. ; . - J « • » 

r y tl  • J *r  m 


•..*  * *h  ir lie-  t»»s  th'-t  the  scuffle iert  of 

T"  T*~”  i 

•"/"i  t.fvf'i  (p  f'=r  ( . ',rid  re  iuc«  it  "ur'/r-r  '.<»  t',-  "onr. 

,.3,.  .>«rL‘  .P-i,  , • r-k  JUtU 

(33,*  fpr(j>c-jrl)  (k  J.)  ( — C ^ / 


.<■;  v>  " -,-r  ”>•  C,  h r>  r.  util  1? 


■>.  « r* 


*.  h .-'  1 ■» 


i/;  j ■ 


whioi  ly  tnak  3f  use  of  h e notation  In*,  "educe?  in  (t*4)  tec  e vritt^n  r 


i8»; 


(-U  r 'Ufi,  n-J  * !)(«,• 


j 

A?  sr.en  Honed  »'*»rore,  when  r l?  even»  the  ri:tti  u'-nd  -hie  «•**  © n :*u: 
reduces  „o  ?sro  end  cancelling  out  the  non-roro  ccrH^ir if-nt  ( 1,!  xr  . »fi,n  jtl) 
frcti  tie  left  hand  side  ns  ©r^resned  1/  (d4)  -re  o.  tain  c,  a**  -t  root  of'  tne  sa.w 
©qua  t .on  (60)  obtained  earlier  by  a different  rr.et.od. 


file  rif’ht  hand  sloe  o tin  equation  (*f.l) , e;.c*?,»t  for  the  f‘.< ©tor 
x 

can  be  written  as 

os)  T (J)(-c,)r"kf j y1  ay, 

flo  5 


s-0 


and  by  tanking  use  of  the  relation 

(3f.)  0-1) 

k-0  T 


it  can  be  reduced  to 


;>1 


n-.i 


(87)  B (r,  jfs+l)c.rfJ*£ 

6 J 


srU 


Using  (84)  and  (85)  we  can,  thus,  write  the  9(.uation  (81)  as 

(88)  B(J+1,  n-j+l)(cj-ig)(r”1)Jfs1&-(-i)r  ^(_i)s(n-J)D(r,  j+8+l)cJr+J'f£ 

L _j  3aQ 


This  eauation  is  to  be  solved  for  c4  to  ret  a minima*  invariant  jrocedure  for 

j 

estimating  F when  the  loss  function  is  given  by  (1).  When  r is  even,  the  factor 
l-(-l)r  — 0 and  we  get  an  equation  of  degree  (r  1).  When  r is  odd,  the  factor 
l-(-l)r  s2  and  the  equation  reduces  to 


(89) 


n-j 

^(-DS(nsJ)B(r,  jfs+Dc  r+j+s  -$B( j+i,  n-j+l)(cr  ii|) 
s=0 


jil.(r-l)*_ 


which  is  an  equation  of  degree  (nir).  In  either  case  there  is  one  and  only  one 
real  root  which  lies  between  0 and  1 and  the  set  of  such  roots  for  j*0, 
minimizes  R- 


An  alternative  way  of  expressing  the  right  hand  side  of  (81)  is  to  rewrite 
(87)  in  the  following  forms 


( 


Q 1 Q 

-l)1*-1  rl  ft(-l)fl(V)  

sstf  ( 1fS+l)(i-4 


j - 


(90) 


= (-l)r~1  rl  y1  (-D8(°-J 

3*0 


Uf  St1)  ( j-Js+2) . . . ( jt-s+r) 

./'j  ,*r  ?2 

• • A *iJfS  ^ ■ -4% 

0 -J0  *0 


fi  (z2 

=(  -1)  rl  \ . . A 3-f  (i-%.  * ^dZv  . . d 4 

Jo  Jo  Jo 

The  equation  (88)  can,  therefore,  also  be  expressed  as 


22 


(91)  B(  j+1,  n-j-ll) (c  - iii)(r_l)  ' 

*>  n-,2 


,l-(-l)rJ(r-l) 


r j r*  ,^2 

. A 2^(1-^)°-^^.  . d* 

0 J0  /□ 


'} 


(92) 


From  this  fora,  it  is  easily  seen  that  for  r«l,  the  e uation  reduces  to 

)o 


B(j+1»  n-j+l)*2\  dZ 


which  shows  that  c,  is  the  red is a of  the  Beta  distribution  with  density 
J 


g(»)=_,T_i — 2j(l^)n_j  , 06H61 

B.j+1,  n-j-tD 


for  J = 0,  1, . . , n, 


as  obtained  earlier  in  Section  4. 


I would  like  to  aekiowledre  sor.e  very  helpful  discussions  with  ’rofessora 
Z.  \!  PirnbsuE  and  H Rabin  durinr  the  preparation  of  this  paper. 
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